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Aldous has introduced a notion of scale-invariant random spatial 
network (SIRSN) as a mathematical formalization of road networks. 
Intuitively, those are random processes that assign a route between 
each pair of points in Euclidean space, while being invariant under 
rotation, translation, and change of scale, and such that the routes 
are not too long and mainly lie on “main roads”. 

The only known example was somewhat artificial since invariance 
had to be added using randomization at the end of the construction. 
We prove that the network of geodesics in the random metric space 
generated by a Poisson line process marked by speeds according to a 
power law is a SIRSN, in any dimension. 

Along the way, we establish bounds comparing Euclidean balls and 
balls for the random metric space. We also prove that in dimension 
more than two, the geodesics have “many directions” near each point 
where they are not straight. 


1. Introduction. Scale-invariant random spatial networks (SIRSNs) are 
a class of random networks defined as a route between each pair of points, 
with three types of properties. First, invariance properties, second, guaran¬ 
tees on mean lengths of routes—in the Euclidean metric—and third guar¬ 
antees on the mean length of intersection of a suitably truncated version of 
the network with a given compact set. It turns out that these conditions 
are enough to deeply constrain the network. For example, all SIRSNs have 
singly-infinite paths for which any subset is included in a route, but no such 
paths are doubly-infinite. 

The only known example is the binary hierarchy model, in two dimensions. 
It consists of minimum-time paths on a dyadic grid where speed depends on 
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the two-valuation. Invariance is obtained by a randomization construction. 
The latter feature is somewhat unsatisfying: the model itself is invariant, 
but realisations exhibit long-range dependence: observation of a small region 
gives much information on the network everywhere. 

A more “natural” candidate for a SIRSN is therefore the Poisson line 
process model. Intuitively, lines are thrown uniformly at random in and 
marked with random speed limits. Slower lines are dense in W^. Then the 
route between two points is the minimum-time path made of segments of 
these lines. Remarkably, even in dimension d>3, when random lines almost 
surely do not intersect, such paths exist. The whole construction is invariant. 
The aim of this paper is to show that the Poisson line process is indeed a 
SIRSN for all d > 2. 

Historically, Aldous (2014) introduced the notion of a SIRSN, and proved 
a number of their properties, including those mentioned in the first para¬ 
graphs of this Introduction. Aldous and Ganesan (2013) give a verbal de¬ 
scription. The motivation was twofold. 

First, Aldous and Kendall (2008) had proved that it was possible to build 
a road network connecting a prescribed set of points that both had routes 
almost as short as the segments between each pair of points, that is, the 
corresponding Euclidean geodesics, and total road length almost as short as 
the Steiner tree, that is, the shortest possible connecting network. However, 
the network was less efficient at small scales. Thanks to their invariance 
properties, SIRSNs have the same efficiency at all scales. It turns out that 
there is a trade-off: for a SIRSN, there is a lower bound on the expected 
total length of the network, which is decreasing in the expected route length 
between two points. 

The second motivation was to give a mathematical abstraction of road net¬ 
works and maps, in particular online maps as they are used today. Namely, 
we may change viewpoint, zoom in or out, and the appearance changes little, 
as smaller roads are shown and hidden. SIRSNs are then defined as statisti¬ 
cally invariant under translation, rotation and change of scale. Moreover, we 
are less interested in the roads than in the routes: how do we drive from A 
to B1 SIRSNs are thus defined by giving routes only, namely unique routes 
connecting pairs of points. It turns out that a notion of “main roads” at 
any scale can be built from the network of routes itself. To wit, the network 
of main roads at scale r would be the network of routes deprived of balls 
of radius r around their endpoints. It is finite in every compact. Similarly, 
Kalapala et al. (2006) have shown that a number of statistics of real road 
networks do not depend on scale. 

Aldous (2014)proved that the binary hierarchy model was a SIRSN, and 
suggested two other possible models for SIRSN, one of which is the Poisson 
line process model. Kendall (2015) has then proved important properties 
of the Poisson line process with appropriate speeds: it does yield a random 
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metric space and this space is a geodesic space. Moreover, in dimension two, 
the geodesics are almost everywhere unique, the geodesics are locally of finite 
mean-length, and the subnetwork obtained from the routes connecting points 
of an independent Poisson point process has finite length in a compact set. 
The latter properties establish a “pre-SIRSN” result, but fall short of the 
full definition. 

As a candidate for a SIRSN, Poisson line process model then fall in a large 
category: networks derived from geodesic spaces. Indeed, from any geodesic 
space, we may build a spatial network by associating to any two points the 
geodesic(s) connecting them. It is not obvious how one might determine 
when such a network is a SIRSN. 

In Section 2, we give a precise definition of a SIRSN and of Poisson line 
processes. We also present other notation and definitions, and recall some 
known results, in particular that the Poisson line process with speed limits 
yields a random metric on W^. We then give a few basic properties of 11- 
paths, that is paths in this metric space. In Section 3, we give sharp bounds 
on the random diameter for this metric of a Euclidean ball, with a few gen¬ 
eralizations. These estimates will be an important tool in several subsequent 
proofs. In Section 4, we prove that geodesics between a given pair of points 
are almost surely unique, in any dimension. Lemma 4.6 will play a central 
role: we introduce the notion of “many directions”, and the lemma states 
that geodesics have many directions at relevant points. A consequence is 
that any path using the same lines as a geodesic will contain these points. 
We will then conclude by noticing that geodesics between the same pair of 
points almost surely use the same lines (Lemma 4.7). In Section 5, we prove 
that geodesics have finite mean Euclidean length. Alternatively, we may see 
the result as supplying a stochastic control of the Euclidean diameter of balls 
for the metric generated by the Poisson line process. Section 6 contains the 
last and most important component of the proof that the Poisson line pro¬ 
cess generates a SIRSN. Intuitively, we establish a sharp control of the total 
length of the intersection of all infinitely many geodesics minus a ball around 
each endpoint, with a given ball. This corresponds to the fact that all these 
geodesics coalesce before hitting the ball and split after leaving it. Bounds 
are given using the pigeon-hole principle and the fact that relevant geodesics 
must use the few fast lines that hit the ball. Einally, Section 7 suggests a 
few potential directions of future inquiry. 

2. Notation, definitions, basic properties. We follow Kendall’s (2015) 
notation whenever possible. 

We write B(x,r) for the ball with center x and radius r. 

We first give the precise definition of a SIRSN. Suppose that A is an 
atom-free measure on a measurable space (A, B). Recall that a Poisson point 
process of intensity measure A is a random set of points such that for any 
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B gB, the number of points N{B) in i? is a Poisson variable with intensity 
A{B), and the number of points N{Bi) in disjoint Bi are independent. Then 
a SIRSN is defined as follows. 

Definition 2.1. A SIRSN is a process that associates to any two points 
xi and X 2 in random routes such that: 

1. Between two specified points xi and X 2 -, there is almost surely a unique 
route TZ{xi,X 2 ) = TZ{x 2 -,xi). It is a finite-length path connecting xi and X 2 - 

2. For a finite number of points xi,...,Xk in consider the subnetwork 
AA(xi ,... ,Xk) formed by the random routes connecting all Xi and Xj. Then 
AA(xi,..., Xk) is statistically invariant under translation, rotation and change 
of scale: if 91 is a Euclidean similarity of then A^(91(3:i),... ,91(xfc)) has 
the same distribution as 9iAf {xi, ..., Xk)- 

3. Let Di be the length of a route between two points at unit distance. 
Then E[Di] < oo. 

4. Let n S N*} be a collection of Poisson processes on with inten¬ 
sity n times Lebesgue, all independent from the SIRSN. Suppose they are 
coupled so that En C En+i- Write H = (JneN* Then the intensity (mean 
length per unit area) p(l) of the following long-distance network is finite: 

IJ (7^(xl,X2)) \ (B(xi,l) UR(x2,1)). 

xi,X2£^ 


Note that Kendall (2015) uses more conditions in his definition, but the 
missing properties are implied by property 4. They were useful to define 
weaker variants of a SIRSN. 

The use of Poisson processes in property 4 makes it look slightly complex, 
but this is a technical shortcut: it allows us to study the network through 
only countably many routes. Morally, we would like property 4 to hold true 
if the long-distance network was defined as the union of all routes between 
all pairs of points of M'^, minus the balls around the endpoints. But there 
would be uncountably many routes, and it would be harder to work with. 

We now turn to the definition of the improper Poisson line process. We 
first need a measure on lines. More details on this kind of process may be 
found in the book by Stoyan, Kendall and Mecke (1996). 

Let be the space of all lines of A line is “un-sensed”, that is, it is 
seen as a subset of without a preferred direction. For K, a compact of 
the hitting set of K is the set of lines that intersect K, denoted as 

[K] = {I e : I hits K}. 

We also denote by the Hausdorff measure of dimension d. With this 
notation, we have the following. 
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Definition 2.2. The invariant line measure ^id{l) is the unique measure 
on the space of lines of that is invariant under Euclidean isometries, 
and normalized by the following requirement: for a compact set K in 
of nonempty interior, the //^-measure of the hitting set of K is half the 
Hausdorff [d— l)-dimensional measure of the boundary of K: 

l^d{[K]) = \md-i{dK). 


The reason for the normalizing constant ^ is to ensure that the measure 
of the hitting set of a flat hypersurface A is its hyperarea md-i{A). 

We will often need the hitting set of a ball, so write ujd-i for the hyperarea 
of the unit sphere Thus ^d{[B{x,r)]) = _ 

Consider the following parameterization of a line h it is given by a direc¬ 
tion and a localization. The direction is an element of the projective space 
It then defines a hyperplane normal to this direction, through a spe¬ 
cial point—the origin—of The localization is a point on this hyperplane. 

The line I is then the line through this point normal to this hyperplane. 

Hence, writing % = for an hyperplane of M'^, we may parameterize 
the set of lines by ®T-L. Notice that if we want to keep track of the 

topology of the set of lines, the product should be twisted, but we only need 
measure-theoretical properties, so we stick to the simpler direct product. 

Recall that the projective space may be seen as the sphere 

with opposite—antipodal—points identified. The projective space then in¬ 
herits the natural metric on the sphere, namely the distance between two 
pairs of antipodal points is the angle between the pair of segments joining 
them, in radians. Up to a null-measure set, the projective space may be more 
simply viewed as a hemisphere. 

With this parameterization, and writing for the unit ball in we 
may write fid as a product measure on 


Mrf = 


1 


md-i®md-i. 


To make a clearer reference to it, we write for 

the probability measure on For a set of lines C, we write 

for the measure of the set of directions of lines in C that go through the 
origin. In particular, for A a subset of and x a point of the value 
“ *]) proportion of directions that appear in the set of 

directions of lines through A and x. We call it the solid angle of A from 
X. Notice that it is not quite the usual definition since pairs of antipodal 
points contribute only once, and since we have normalized to a probability 
measure. This is because we are concerned with orientations of undirected 
lines, rather than angles of directed lines. 
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A special case we shall need is the measure of a cone, that is, all the 
lines with angle less than 0q from a given line through the origin. This 
measure is proportional to the area on the hemisphere hit by the lines, 
hence proportional to sin'^ ^ 9 cos 9d0 = sin'^ ^(9o)/(d—l). The integral 
up to 9o = 7r/ 2 has value 1, so that 

( 1 ) /r^'^j^(cone of aperture Bg) = sm^~^(9o). 

The uniform Poisson line process is the image of the Poisson point process 
on with intensity measure 

We define our improper Poisson line process by adding a mark on each 
line, a speed limit. Namely, the improper Poisson line process is the image 
of a Poisson point process on PW^~^ ®T-L ( 8 )M+ with intensity measure 
given by the density 

dyLda{^,v) = dfj.d{l){7 - du 

for 7 > d. Kendall (2015) does define this process for all 7 > 1, but the 
relevant case for SIRSNs is that of 7 > d. 

In words, we have more and more slower lines, following a power law. 
Since f v~'^ du diverges at zero, the lines are dense in However, lines 
faster than any given speed are not dense. In particular, the number of lines 
faster than vq hitting a convex set K is a Poisson variable with parameter 

roo 

fJ-d,'r v) :l£[K] and v> Vo) = UdiiK]) / {-f - l)v~^ dv 

J VQ 

= ^md-iidK)vQ^'^~^K 

We call n = n(d, 7 ) the corresponding random process of marked lines 
{l,v). Since the dimension d and parameter 7 will always be clear from 
context, we drop them in the notation. Notice that the total number of lines 
is almost surely countable. If {l,v) G H, we say that the speed of line I is v 
and denote it v{l). 

For a subset £ of lines, we write H^ for the restriction of H to these lines, 
that is. Tic = ^ III particular, the line process restricted to lines 

hitting A but not B is 

We denote S the silhouette of H, that is the random set in made of 
all the lines of H, that is 5 = {x G : 3{l, u) G H : x G /}. We also write 5^^ 
for the random closed set in made of all the lines {l,v) in H such that 

V >Vo. 

We may then define H-paths. 

Definition 2.3. A finite-time H-path is a locally Lipschitz path in 
respecting the speed limits imposed by H. More precisely, it is a continuous 
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M'^-valued function 

C:[0,T]^R^, 

with T finite, such that for almost all t G [0,T], either: 

• the speed is zero: ^'{t) = 0; 

• or the path follows a line in IT: there is a u > |C^(t)| such that {^{t) + 

gH. 

We call T = T{^) the time length of the path or just its time for short. 
An infinite-time Il-path ^ is the same, with T(^) = oo, except that its 
domain is [0, oo). 

A Il-path is a finite or infinite-time Il-path. 

In an abuse of notation, we write ^ G IT. 

Notice that the image of a Il-path is not necessarily contained in the 
silhouette S. It only needs to have speed zero outside S. The remark is 
especially relevant in dimension at least three, where the lines never cross. 
However, since the lines are dense in it turns out that there are paths 
joining any pair of points, without any segment in \ 5. We give a clearer 
intuition of their tree-like structure in the proof of Theorem 3.1 and Figure 1. 

We write H"^ = {^ G H : ^(0) = a and = b} for any two points a 

and b. 

The following theorem is a union of results from Kendall’s (2015) paper. 

Theorem 2.4. Almost surely, all finite-time H-paths have finite Eu¬ 
clidean length. 

Almost surely, there are finite-time H-paths between each pair of points of 
M'^. Moreover, for any two points a and b, the infimum Tab of time lengths 
T{ff) ofH-paths f G is attained. 

Hence, with the metric d{a, b) = Tab is a random metric geodesic space. 

We call this metric time length or H-length. Time diameters and similar 
notions are defined in the same way. 

We denote by Af the random network made by all the geodesics connecting 
all pairs of points in Our aim is to show that W is a SIRSN. 

It is often possible to define similar metrics on other sets T of marked 
lines, though they might not be geodesic. We then speak of T-length. The 
typical case is when we restrict H to a subset of lines C, yielding T = H^ 
and n^-length. 

We now introduce some notation and remarks to make easier manipulat¬ 
ing paths and geodesics: 

• If^G^“^ we often write f,ab instead. 


J. KAHN 








Fig. 1. We cover each hall Bi with a ri-net, and connect each pair of balls of the net 
with a segment, he it a point (a). We then do the same at each scale, for rn+i-nets of 
the balls of the Vn-nets (b). Each point belongs to a ball of the ri-net, so we build a path 
between two points starting from the segment connecting their balls (c). We then (d) and 
(e) connect the points to the segment endpoints with the segment connecting their balls in 
their respective r 2 -nets. At stage n, the path (f) is made of 2" segments at each scale n. 
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• Concatenation of Il-paths is denoted ^ac = CabCbc, that is, ^ac(t) = Cab(t) if 
t < T{^ab), and ^ac{t) = ^bc{t - T{^ab)) if t > T{^ab)- 

• We use the letter g for geodesics, and usually gab for a geodesic from a to 
b. We say that gab is unique if there is a unique geodesic from a to b. 

• In an abuse of notation, we identify a Il-path ^ with its image in 
whenever it is clear. Hence, we may write x £ ^ if there is a time t such 
that ^{t) =x. Similarly, if gab is unique and c,d£ gab, we say that gcd is 
included in gab- 

• If gab is unique and ce gab, then gab = gacQcb- 

• T {gab) = Tab- 

• For a line I £ H and a H-path ^ : [0,r] —>• we define the intersection 

length of ^ and I as L^{1) = mi{l n ^([0,T])). 

• If the intersection length of ^ and I is not zero, we say that I is in the 
support of or that it is supporting We write I C Moreover, we 
denote the support of ^ by = {/ £ H : Lg(/) > 0}. 

• We define the intersection time of I and ^ as T^{1) = mi{[0,T] n^“^(/)). 

• In particular, almost surely, for all geodesics g, we have the following 
equality and decompositions: 

Lg{l)=v{l)Tg{l) foralU£n, 

(2) T{g)=Y,Tg{l) 

leCg 

= J2v{l)Lg{l)- 

l&Cg 

• If the support of ^ is included in /I, that is C C, we say that ^ is 
n^-path. We abuse notation by writing ^ £ H^. 

• Similarly, we write ^ £ H^f if ^ £ H^ and ^ . 

Finally, a few more generic notation. We call internal e-net of a subset A 
of any maximal subset xi,. -. ,Xk oi A such that \xi — Xj \ >e for all i^ j- 
For a set A C M'^, we write A^ for its r-widening, that is the Minkowski 
sum H'’ = H0 H(0,r). We denote the maximal speed of a set of lines C 
by Umax(^)- That is, Umax(T) = supi^£^v{l). Notice that if those lines all 
hit a compact set K, this supremum is a maximum. We abuse notation 
by writing Umax(^) = ^^max([^]) for AeR^. We use C, c, ci,..., for any 
positive constant. 

3. Il-diameters of sets. We start with giving a more quantitative version 
of Theorem 3.6 in Kendall’s (2015) article. Namely, we show that in a given 
precompact set A, any two points can be joined in finite time, and that the 
largest time between two such points—the H-diameter of A —is not too big 
with high probability: this random variable has more than an exponential 
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moment. We include a generalization that we will need later on, by allowing 
the possibility of ignoring lines hitting forbidden areas F. 

Theorem 3.1. Recall that d>2 and 7 > d. Let amin = . Let 

LI < rimax = o-^d define Omax = I A. Choose a sueh that 

ctmin < a < a m ax . and note in partieular that a > 1. 

Let A and F be two subsets sueh that, for some r > 0: 

• A is eonnected. 

• From any x € the solid angle hl^li{[F — x]) of F is less than LI. 

o 

• A may be eovered by N balls Bi of radius r. 

Then there is a Ti depending only on a, LI, ^ and d such that, for any 

o , 

Emax = l/(2A^(2a + 1)“) > e > 0, with probability 1 — e/emax; the diameter of 
A on nj^i./(c.-i)]y[p.] satisfies the bound 

Ta,f = sup inf r(0 

(3) <r^^(^ln-j with 

T o = 

r,N 

In particular, this maximal time has all exponential moments, and more: for 
any 6 <T y*'"’’ we have 

r,N 

(4) E[exp(5r^^y] < 00 . 

The proof is a slight variation on that of Theorem 3.6 in Kendall’s (2015) 
article. 

o 

Proof of Theorem 3.1. Since A is connected and covered by N open 
balls Bi of radius r, we may build a path between any two points of A by 

o 

concatenating at most N paths between two points of A belonging to the 
same ball Bi of the cover. 

We now recursively build a path between each pair of points xq and yo of 
Bi, in a binary tree-like fashion. First, we specify 

rn = ra~'^. 

We will choose corresponding speed limits later. Given such Vn- 

• We call Aq the set of balls 
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• To any ball G An, we associate an internal r„+i-net of that ball. It 

may be viewed as a collection of balls of radius r„+i. We then 

define An+i as the set of all these balls for all together. 

• We have thus built nested internal r„-nets. 

• For any two balls and belonging to the internal r„+i-net 

of the same ball B^^') in An, we find a line of speed at least Vn that hits 
both and Bj^~^^\ but not F. We will have to prove this is possible 

with high enough probability. 

• For Xn and both belonging to B^'^^ G An, we may then hnd two points 

Xn+i and Un+i such that: and Xn+i (resp., and yn+i) belong to the 

same ball G An (resp., and Xn+i and t/n+i both belong to 

the same line of speed at least Vn- 

• We may then build a path between Xn and i/n as a concatenation of three 
paths: CxnVn = ^xnXn+i^xn+iyn+iCyn+iyn- The middle one is a segment. The 
other two are paths between points of the same ball in An+i- 

As illustrated in Figure 1, we thus obtain a path between xq and yo that 
is made of exactly 2 ” segments for each n, each at speed at least Vn between 
two balls of the same internal r„+i-net of a ball of radius r„. Moreover, since 
the points xq and yo are in A, all of the segments are between points of the 
Minkowski sum i?(0, r^) = 

This construction has built a path for each pair of points xq and yo in 
Bi. Since segments between balls of the same internal r„_|_i-net of a ball of 
radius are at most ( 2 r„ + 2 r„+i) long, and are at speed at least Vn, the 
nj^r/(a-i)]y[p']-diameter Ta,f of A is bounded from above by 

(5) 

n=l 

on the event that this construction is possible. 

Now to control the probability of this event, we need: 

• a bound on the number of pairs of balls in the same r„-net; 

• a bound on the probability that the fastest line hitting two such balls but 
not F is slower than Vn- 

We obtain the first bound by using the formula Indeed, 

each r„-net of a r„_i-ball is then the same as a a“^-net of a radius 1 ball. 
Since the balls B{si,l/{2a)) centered on the points of such a a“^-net are 
disjoint, and all included in a ball of radius 1 + l/( 2 a), a volume argument 
shows that there at most (2q; + balls in each r„-net. So that there are at 
most N{2a + 1)'^"’ balls in An, and at most N{2a + pairs of balls in 

the same r„-net. 
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Fig. 2. Lines that hit both balls B{x,rn+i) and B{y,rn+i) but not F (not necessarily 
connected) contain lines through z £ D{x,rn+-i_l2) and B{z + x — y,rn+i/2). These are 
included in the cone of angle 9o through z. 

We now consider these lines that hit both of two r„+i-balls B{x,rn+i) 
and B{y,rn+i) in an internal r^+i-net of a r„-ball, but that do not hit F. 
We have seen that they were in We may then use the hypothesis 

of the theorem on the solid angle of F. 

We want a bound on 

Fd{[B{x,rn+i)] n [B{y,rn+i)] n [F]'"). 

Now, as illustrated in Figure 2, let D{x, rn+i/2) be the hyperdisk centered 
on X with radius r„+i/2 and included in the hyperplane perpendicular to the 
vector x — y.lt is included in B{x,rn+i)- Moreover, for any z G D{x,rn+i/2), 
the ball i?(z + y — x,r„+i/2) is included in i?(y,r„+i). 

Now we may write 

Fdi[B{x,rn+i)] n [B{y,rn+i)] n [F]'^) 

> md_i(F(x,r„+i/ 2 )) 

X ^.inf y^l^{[B{v-x,rn+i/2)]f\[F - zf) 
z^D{x,rn+il^) 

>md-i{D{x,rn+i/2)) x {y'°l^{[B{y - x,rn+i/2)]) - Q). 

Now, y — X < 2r„ = 2ar„+i. So that the lines in [B{y — x,r„+i/2)] are in 
a cone of angle Oq such that sin(0o) < (4a)“^. So that by formula (1): 

(6) M2i(|B{</-i,>'„+i/2)|)<(4a)‘-‘'. 

Since a < Omax) we have > 0. 

Moreover, with Vd-i the {d — l)-volume of a unit {d — l)-ball, we have 

/ \ d—1 

md-i{D{x,rn+i/2)) = > 

/r,([F(x,r„+i)]n[F(y,r„+i)]n[F]^) < - 0) ' 
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using the notation 0^,0 to summarize and emphasize that this bound de¬ 
pends on a and but not on n, r, or any other feature of A and F. 

We thus get 

/ \ 

IPbmax([-B(a;,rn+i)] n [B{y,rn+i)] n [F]'^) < Vn] < exp|^-Ca,n-^J. 

Multiplying by the number of relevant pairs of ball, we then obtain that 
the construction is possible except on an event of probability at most 


( 7 ) 


^iV(2a + l)'^(^+i)exp 


n=0 


—Can 


rd-i 
' n+1 
7—1 

vA 


Taking 


Vr,. = 


(d-l)/(7-l) 

n-l-1 


((n + l)c„,f,ln(l/e))V(7-i) 
and replacing in bound (7), we see that this becomes less than 

00 

^ N{2a + l)d{n+i)^n+i < 


for e < l/{2N{2a + l)'='(^+i)) = so that ffmax does not depend on A or 

o 

F except through a and N. 

Replacing Vn and with their value in the bound (5), we get 


Ta,f < iV^2”+^(l + 


a 


Jn+l 


n=l 


= iV ( (1 + 


X ( In 


n=2 

l\ \ l/(7-l) 


where the hrst parentheses correspond to Ti and Ti is hnite since a> a„ 
We have thus proved formula (3) of the theorem. 

The moment (4) is simple integration: 


1 F£max / 

E[exp((5r^”p^)] <C-\ -/ exp( S 

’ ^max Jo \ 


<c + 


Ff 

£max Jo 


r,N 

Emax 


]^\ l/(7-l)\ 7-1 

In- I 

c 


de 


r,N 


de 


< 00. 
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Here, we use 5 <T and C corresponds to the integral between e. 


r,N 


max 


and 1 , bounded by the value of Ta^f for e = ^max- D 


Remark 3.1. If no lines are forbidden, that is if R = 0 , then a can be 
taken as big as we wish, so that the lines used all hit as small a widening of 
A as we want. 

We may slightly generalize the theorem: instead of using a forbidden area 
F, we could use different conditions for which lines to accept. The important 
property is that we must have enough relevant lines hitting pairs of balls in 
a r^+i-net. 

There are a few optimisations that could be used to gain slightly in the 
constants. For example, we have written the proof with one subnet inside 
each ball. If we had used a single r„-net of a correctly widened A, we would 
have only about pairs of ball, allowing a bigger ^max- The result 

stays essentially the same, however, as proven by the following proposition. 

Proposition 3.2. For any two points x and y, their U-distance T^y 
does not have a moment with higher exponent on the time Txy than in for¬ 
mula (4)- That is, for any 5 > 0, for any rj> ^ — 1, we have 


]E[exp(5T,^j^)] = oo. 


Proof. Say that \x — y\=r. A path from x to y has to go from x to 
the border of B{x,r). Hence, it must use lines hitting B{x,r) for a length 
at least r. So that T^y is controlled by the fastest line hitting B{x,r). Now, 


Pbmax(H(x, r)) < Ve] = exp(-c(r)u^ '^) = ^ 


with 



where c(r) depends only on r. 

So that on an event of probability at least e, we have the bound Txy > 
r/ve = C 2 (r)(ln(l/e))^/('>'“^\ Hence, for some positive constant C 3 depending 
on 6, rj and r: 



= 00 


since r //(7 — 1 ) > 1 . □ 
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4. Almost sure uniqueness of Il-geodesics. In this section, we prove that 
the network Af satisfies the property 1 of a SIRSN, that is that between two 
specified points x and y in there is almost surely only a unique route. 
Since routes are the Il-geodesics, this is equivalent to almost sure uniqueness 
of the geodesic gxy 

The case in dimension 2 has already been established in Kendall’s (2015) 
paper. The following proof, on the other hand works in all dimensions more 
than 2, as stated in Theorem 4.8. Note that it does not work in dimension 

2 . 

The strategy is the following: 

• We introduce a concept of many directions, with the following property. If 
a Il-path has many directions near a point x, any finite Il-path supported 
by the same lines contains x. 

• We show that almost surely all geodesics have many directions near all 
the ends of their component segments, except for the two extremal points. 
This is the step where d > 3 is needed. 

• We show that for specihed x and y in almost surely all geodesics from 
X to y are supported exactly on the same lines. 

• So that almost surely all such geodesics contain the same segment ends, 
and this will prove they are the same. 

The author thinks the proof is very technical for something that looks 
clear enough, but could not find an easier way. Maybe the need to work 
with tree-like paths in dimension at least three is the reason why there is no 
obvious argument. Hopefully, the concept of many directions can be useful 
elsewhere. 

We first state two technical lemmas we need for the proofs in Section 4.1. 
We introduce the notion of “many directions” and give some cases where 
paths have many directions in Section 4.2, culminating in Lemma 4.6. We 
prove that geodesics must use the same lines in Section 4.3, and end the 
proof of uniqueness in Section 4.4. 


4.1. Technical lemmas. The first lemma yields a control on the propor¬ 
tion of balls in a nested set that are hit by lines faster than a threshold 
appropriately scaling with their size. 


Lemma 4.1. Let d>2. Let a > 1 be a scale factor. We write p = < 

1. Let Bi = B{xi,ri) for 1 <i <n be a set of nested halls with ri = rQlW for 
some tq. For some vq, define 


Vi = , 


Vo = {{l,v) G n : u > Uo}, 

Vi = { il , v ) G n : Vi-i > v > Vi }. 
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Define now Hi as the number of balls Bi^k smaller than Bi that are hit by 
lines faster than Vi, and In as the number of balls Bi that are hit by lines 
faster than Vi-i, that is, 

Hi = if{k > 0 : Vj n [-Bj+fc] 7 ^ 0 }, 

n—1 

In = #{* < n : 3A: > 0 : [Bi\ n Vi-k 7 ^ 0} < ^ Hi. 

1=0 

Then A = E[ 7 ^{Vj Pi [i?*]}] = (1 — p)u}d-irQ~^ / {2 vq~^) independently of i>l. 
Then, for any 5 > 0; 

(8) P[/n >6n]<(^{l + \5) V) "• 

Precise but more cumbersome bounds are given in equations (11) and (12). 


Proof. We have 

E[#{Vi n [B,+k]}] = E[#{V* n [Si]}]/ 

= Xp^ for all i > 1 and k>0, 

00 

E[#{Von[Sfc]}] = A/ 

m=0 


1-p 


Xp^ for all A: > 0. 


Thus, 


F[Hi = 0 ] < 1 for all i > 0 , 

F[Hi = k>0]= P[Vi n [Sj+fe] 7 ^ 0 and Vi n [Sj+fc+i] = 0] 

= exp(-A/+^) - exp(-A/) using [Sj+i] C [Bi] 

< l-exp(-A(/-p^+i)) 

< Ap^(l — p) for i>l. 

Similarly, 

Apfc+i \ / \pk 

- exp - 


P[So = fe > 0] = exp ( — 
< A/ 


l-p 


V 1 -p 


Hence, for all 0 < a < 1/p: 

00 

E[a^'] < l + ^A/(l-p)a^ 

k=l 
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(9) 


( 10 ) 


= 1 + A(1 -p) 


ap 


1 — ap 


for all i>l, 


E[a^o] < 1 + ^ Ap' 




k=l 


= 1 + A: 


ap 


' 1 ) 

1 — ap 

E[a'^"] < E[a^i‘=o^ V 1 
< fl + A(l-p): 


ap 


1 + ap(A - 1) 


1 — ap J 1 + ap(A(l — p) — 1) 


< 1 + A 


ap 


1 — ap 


We now use the Markov inequality: 

( 11 ) 

and optimize upon a. Exact optimization requires solving a degree-two equa¬ 
tion and yields a cumbersome solution, so we shall only use here the solution 
for the second approximation when A = 1, that is a = 5/{{l + S)p), so that 
ap/(l — ap) = 6. Using (10), we get 


P[/„ > 5n] < 


1 + Xap/ (1 — ap) 

1 -L A(1 -p)ap/(l - op) 


l-LA(l-p) 


ap 


1 — op 


^—Sn 


( 12 ) 


1 + A<5 


1 -L A(1 -p)5 


(l + A(l-p)5) 


1 + 5 


p 


< 


(1 + A5) 



□ 


The second lemma gives a guarantee that except on exceedingly rare 
events; no significant fraction of uniformly randomly oriented lines are clumped 
together in a small number of directions. 


Lemma 4.2. Let {X,d,fj,) be a metric space with a probability measure 
p, and such that for all positive and small enough s, for any point x in X, 
the volume of the ball is bounded in this way: 

(13) cie'^~^ < p[B{x,e)] < C 2 e'^“^, 

for Cl and C 2 constants depending only on the space {X,d,p). 

Let a > /3 > 5 > 0. Consider {sjjjg/ a n~'^-net of X. Consider {xj}j<an 
an random p-i.i.d. points in X. 
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Then there is no subset of j3n points in {xj} that are all contained in at 
most 6 n balls of the net B{si,n~'^), except on an event of sub-exponential 
probability, at most with C depending only on a, {3 

and 6 , ci and C 2 - 

Moreover, the projective space {PW^~^, 9, satisfies the hypotheses 

for {X,d,p). Here 6 {li,l 2 ) is the angle between two lines, and hlili 
natural probability measure on defined in Section 2. 

Proof. Since all the points in the n“^-net are in disjoint balls of radius 
n~'^l2, there are at most points in the net, for n big enough. Now, 

with C 4 ,C 5 ,C 6 depending on a, (3 and 5: 

P[There are f3n points Xj all contained in 5n balls B{si,n~'^) of the net] 

< ^{subsets of fin points} {subsets of 6 n balls} 

X ¥[f3n specific points are all contained in 6 n specific balls] 

In the projective space, where points are seen as lines through the origin 
in M'^, a ball of radius 9q is exactly a cone of aperture Oq. We recall formula (1) 
and 9/2 < sin(0) < 9 for positive small 9, and we conclude that property (13) 
is indeed satisfied. □ 

4.2. Many directions. Having many directions near a point x intuitively 
means that the lines used near the point have so many different unit vectors 
that the only way to touch all those lines (a tour) with a finite curve is by 
touching most of them near x. 

Definition 4.3. For a set of lines C, = {Ijjj^j and a subset of the Eu¬ 
clidean space X C M'^, a C-tour in X is a, curve / in X such that for all 
j G J, there is a tj such that f{tj) G Ij. If / is rectihable, the tour is said to 
be finite; else it is infinite. 

Recalling the notation for the support of f: 

Definition 4.4. A finite H-path f has many directions near a point x 
if, for all e > 0, all /l^-tours in \ B{x,s) are inhnite. 
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As a remark, this concept is only interesting in dimension at least three: 
in dimension two, a circle is usually a tour, and it is finite. 

Proposition 4.5. Let two finite H-paths ^ and rj. If ^ has many direc¬ 
tions near a point x G and its support is included in that of rj, that is, 
C Crj, then x G rj. 

Proof. By Theorem 2.4, the finite Il-path r] has finite Euclidean length. 
Moreover, C Cr], so that r/ is a £^-tour in Since f has many directions 
near x, there is no finite £^-tour in \ B{x,e), for any e > 0. So that 
B{x,e) r\r] 0. Any finite Il-path is closed in hence x G rj. □ 

Lemma 4.6. Letd>3. 

Let I be a prescribed line independent of 11. Almost surely, for all x Gl, 
for any y ^ I, all geodesics Qxy have many directions near x. 

Proof. It is sufficient to show that the claim holds for all x G lriB{0, R), 
for all R 0. It is s^lso sufficiGnt to work with o, fixod £, cinci prove thcit for ciny 
y the e-widening of I, for any geodesic pxy, all -tours on M.‘^\B{x,e) 
have infinite length. 

Let = rofa"^ for ro = e and an a > 1 big enough to be determined later. 
We consider one-dimensional internal ^-nets Nn of B{0,R) Dl. 

In particular, there are at most — + 1 points in Nn- 
Now, for any point in Nn, we may build a set of nested balls {Bi}i<n 
depending on with the following properties: 

• The ball Bn is centred on Xn- 

• All balls Bi = B{xi, Vi) are centred on a point Xi in the net Aj. 

• The centres lie deep inside the previous ball, that is Xi G B{xi-i,^-^). 
This stems from the fact that Aj_i is an ^^^-net. 

We also use the notation of Lemma 4.1 for Vi and V*. We have to choose a 
good speed vq. 

We then show that the nested balls have the following properties, with a 
T to be determined later, except on an event of probability o(r“^): 

1. There are at least ^ balls Bi that no line faster than Uj_i hits. 

2. There are at least ^ indices i such that the time diameter of B{xi,^)\ 

B{xi,'^) on is less than . 

3. There are at least ^ balls Bi that are hit by at most r lines of speed 
between Vi and Uj_i, for some fixed r. Moreover, all those lines have inde¬ 
pendent, uniformly random directions. 
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The T will yield the uq we need to continue the proof, that is, 


(14) 



We deal with these properties in reverse order, and first consider prop¬ 
erty 3. 

Since the sets of lines Vi are disjoint, the events that more than r lines in 
Vi hit Bi are mutually independent (in i), as are their directions. Moreover, 
since the Bi are balls, the direction of the lines of the isotropic Poisson line 
process that hit it are uniformly random. 

So that the number of such events is a binomial random variable 
B{n, q) where q is the probability of a single event. The Chernoff bound for 
a binomial is, for 6 > q, 


F[Xn > 6n] < 




n 


Here, it suffices to take 5 = |. For a big enough, if, for example, we take 
q < a~'^, the above bound is negligible with respect to a~^. 

Now the probability g of a single event is the probability that a Poisson 
random variable with parameter not depending on i is bigger than r. We 
then just choose r so that q < a~^. 

Now for property 2. 

Since the sets of lines [Bi] \ [Hj+i] are disjoint, each of the events that the 
time diameter of B{xi, ^)\B(xi, y) on is less than 

are independent. Arguing as for property 3, we merely have to choose T such 
that the probability g of a single event is less than a~'^. 

Now, by scaling, the sets Ai = B{xi,^)\B{xi,^) may all be covered by 

o 

the same number N of balls with radius y, and centres in A,. Moreover, 
Ai is connected. Now the y-widening of Ai is included in Bi, and since 
Xj+i G i?(xj, y), we know that Xj+i is at distance at least | of this widening. 
So that the maximum solid angle of Hj+i viewed from any point of the 
widening is a decreasing function of a, hitting zero when a goes to inhnity. 
Thus, for a big enough, we may apply Lemma 3.1 with H < Hmax- Replacing 
e in bound (3) with l/(emax«~^) to ensure q < we find that property 2 
is ensured if we choose 


T = NTi 



(7-d)/(7-l) 


(ln(Q:^/emax))^^^^ 


= 0(ln^/^'’' 

This choice of T ensures that vq = 0(ln“^/^'’'“^) a). 
As for property 1, we apply Lemma 4.1. 
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We substitute <5 = | in equation (8), and use p = '^ < a ^—since d > 

3—and A = 0{vq~^) = O(lna). We obtain 

2n 

For less than — indices i, no line faster than Uj_i hits the ball Bi 
5 


< (cln(a)a 
= o{a-^). 


Let us now consider a geodesic gxy from a point y outside Bq to a point 
X inside Bn- Since the balls are nested, they are all crossed by the geodesic. 
Even the sets Ai = B{xi,^)\B{xi,^-) introduced in the proof of property 2 
are all crossed, that is we must pass from the spherical boundary S{xi,^) 
to the smaller boundary iS'(xj, y). 

If the three properties are satisfied, then there are at least indices 
i for which the conditions are simultaneously satished. For such an i, by 
property 2 the geodesic will go from S{xi,^) to S{xi,^) in time at most 
Since the two boundaries are ^ apart, the geodesic must 
use a line with speed more than = voa = Vi within 

Ai- Now, Ai C Bi and, by property 1, there is no line faster than Uj_i in Bi- 
So that the geodesic must use a line in Vi H [Bi ], whose cardinal is less than 
r by property 3. Since the Vi are disjoint, we have proved that the geodesic 
must use at least ^ distinct lines among a set of size at most rn of lines 
with uniformly random direction. 

We may then apply Lemma 4.2. 

We take 1> rj > l/{d — 1) and 6 = 1/20. With probability 1 — 0(n“'^"'), 
all geodesics from a point outside Bq to a point inside Bn must use at least 
n /20 lines (depending on the geodesic), each with a direction in a different 
ball of a n“^-net of the projective space. Since there are at most a constant 
number ci of points in a n“^-net at distance less than 3n~^ from a given 
point, we may choose n/( 20 ci) of those lines, so that each pair of them 
makes an angle at least n~^. 

Let us consider a fixed vq = e. Since d{xi,Xi-i) < and x G B{xn,rn), 
and since r* = = o{en~^) for i> C 2 In re, we know that Bi C B(x, ere“^/3) 

for i > C2lnre. Among our re/(20ci) lines, at most r hit each Bi- Hence, 
up to removing calnre of our re/( 20 ci) lines, all those lines hit the ball 
B{x,en~'^/3). So that, by elementary geometry illustrated in Figure 3, for 
any two lines h and I 2 in our collection, no point of li \ B{x,e) is closer to 
a point of I 2 \ B{x,e) than en~'^ jA. 

We have thus proved that for any geodesic gxy with x G Bn and y ^ Bq, 
any Cg^y-tom in B{x,s) must contain c^n points that are c^n~^ apart 
pairwise. Hence, it has length at least cere^”^, going to infinity with re. 
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Fig. 3. We project on the vector space spanned by h and I 2 , with h parallel to the 
first coordinate and x as the origin. All points of h then have second coordinate at most 
n~'^e/3 in absolute value. On the other hand, let A be a point in I 2 f] B{x,n~'^e/3) and 
B a point in I 2 \ B{x,e). Then the vector A^ has absolute second coordinate at least 
(e — n“''e/3) sin(n“'') > for large n, so that B has absolute second coordinate at 

least . 

All those properties were obtained except on a set of probability o{a~"'). 
Since the net Nn has cardinal of order a^, this is true simultaneously for all 
sets of nested balls built on all the points Xn in Nn, except on an event of 
probability going to zero. 

Letting n go to inhnity yields the lemma. □ 

4.3. Geodesics use the same lines. 

Theorem 4.7. For all x and y in almost surely all geodesics Qxy 
are supported on exactly the same lines. 

Proof. A similar result is used during the proof of Theorem 4.4 in 
Kendall’s (2015) article. The idea of the proof is to “slow down” the lines 
not used by a specific geodesic. Then all geodesics that use one of those 
lines become slower and are no longer geodesics. The set of slowed speeds 
has infinite measure relative to the initial speeds. 

Almost surely, there exists a finite R such that all geodesics from x to y 
are included in B{0,R), by Theorem 2.6 of Kendall. Almost surely, for any 
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i?, there are countably many lines intersecting B{0,R), all with different 
speeds > ^2 > • • • (the information here is that there is a measurable 
order-preserving bijection from N to the inverse line speeds). We then write 
li for the line with speed Uj. 

Let i G N and let us fix all the other speeds Vj for j i. Let us suppose 
there is a speed between Vi-i and Vi+i such that if Vi = v^, then there 
are at least two geodesics and g^y from x to y, one of which is supported 
by li, and the other not. That is k C gl-y and k [Zi g^y. 

Now, the time length of a path not supported by li does not depend on Vi. 
So that g'^y is the fastest of these paths for all Vi, with constant time. On the 
other hand, by decomposition (2), the time length of g\.y is decreasing in Vi. 
Hence, if Vi > , all geodesics from x to y are supported by li. Conversely, 

if Vi < v^, no path supported by k is a geodesic from x to y. Indeed, such a 
path would be as fast as g'^y at speed Vi, hence faster at speed v^. 

So that, when the other speeds are fixed, there is at most one value 
of Vi such that k is in the support of some geodesic from x to y and not in 
the support of another such geodesic. 

Now, we may disintegrate the measures on the line speeds. For almost all 
fixed line speeds for all lines except k, the measure /r* for Vi has a density 
(namely, it is proportional to J. Hence, g-i{{v^}) = 0. So 

that almost surely, either kn g for all geodesics g from x to y, or li\^ g for 
any such geodesic. This is true for all lines li, completing the proof. □ 

4.4. Almost sure uniqueness. We may now state this section’s main re¬ 
sult. 

Theorem 4.8. In dimension d>2, for all x and y in M'’*, almost surely 
there is a unique geodesic gxy from x to y. 

Proof. Kendall (2015) proved the dimension 2 case in Theorem 4.4. 
We then assume from here on that d>3. 

Consider any line / G H. We apply a Palm distribution argument: since 
H is a Poisson process, H \ {?} still has the same distribution as H. So that 
almost surely, by Lemma 4.6, all geodesics on H \ {1} with an endpoint 2 : G I 
have many directions near 2 . Since the lines are countable, this is true for 
all lines simultaneously. 

The rigorous way of writing the former paragraph is through Slivnyak- 
Mecke formula [originally proved by Slivnyak (1962); see, e.g., 
Mpller and Waagepetersen’s (2003) book. Theorem 3.2, for a modern treat¬ 
ment]. A Poisson point process X with intensity g takes value in the set 
of locally finite point configurations Nif. For any nonnegative measurable 
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function h on x Nif, we have 


E 


= [ 

■xex J 


= / E,[h{x, X)] dfi{x). 


We apply the formula with the underlying point process for 11, so that 
X = {l,v) are marked lines. The function h is the indicator function 

{ 1, if 3z G I, y ^ I, gzy Il-geodesic such that 

Qzy does not have many directions near z, 

0 , otherwise. 

Lemma 4.6 then ensnres that the expectation in the integrand in the 
right-hand side is uniformly zero, so that, almost surely, for all I G 11, all 
geodesics on 11 \ {^} with an endpoint z G I have many directions near z. 

Let gxy be a geodesic from x to y. Let {zi} be the set of endpoints of all 
the segments of the geodesic gxy except x and y. For any Zi, since it is a 
segment endpoint (on line 1), there is a tt snch that guzi or gz^u is included 
in gxy, and such that I is not in the support of this sub-geodesic. Hence, this 
sub-geodesic is also a geodesic in n\ {/}, and has many directions near Zi. 
A fortiori, gxy has many directions near Zi. Moreover, by Lemma 4.7 almost 
surely all other geodesics from x to y have the same support. Hence, all 
these geodesics include all the segment endpoints of all the geodesics from 
X to y. 

Now two geodesics from x to y must pass through the same points in 
the same order: indeed if gly = gladlbafy and gly = glbOlaaly, then giggly or 
9xb9by would be shorter than both. 

So that all geodesics from x to y pass through their segment endpoints in 
the same order, so they are the same. □ 


A byproduct of the proof is the following remark. 

Corollary 4.9. Let d>3. Almost surely, for any point x not on a line 
o/H, that is x^S, all geodesics containing x have many directions near x. 

Proof. We use the step in the former proof, that almost surely, for all 
I G H, all geodesics on H \ {^} with an endpoint G I have many directions 
near z. 

For any e, since x ^ S, the geodesic to x will leave a line of H at a point 
2 ; G B{x,ef2). Now the restriction of the geodesic to gzx has many directions 
near z. Since B{z,ej2) C B{x,e), g has also many directions near x. □ 
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5. Geodesic length has finite expectation. This is property 3 of a SIRSN. 

Theorem 5.1. Let £ M'^. Then the Euclidean length of the H-geodesic 
between x and y has finite expectation 

E[Lxy] < oo. 

Proof. For any r, any Il-path containing x whose Euclidean length is 
more than r has to intersect B{x,r) on a length at least r. Hence, it must 
spend time at least r/vmax{B{x,r)) in that ball. If the geodesic gxy has 
Euclidean length Lxy and time length Txy, we thus obtain the constraint 

(15) Vraso,{B{x, r)) > v/Txy for all r < Lxy. 

By equation (3) of Theorem 3.1, there is a T such that with probability 
at least 1 — 

(16) Txy<{n + lf/^^-^'^T = Tn. 

We now consider a collection of radii 

D = 2 Vo 

for I between 0 and m, with tq and m to be chosen later, possibly depending 
on n, but independent of H. 

If Lxy > Vm, then Lxy > ri for all I £ [0, m], so that constraint (15) must be 
satisfied for each r^. In particular, on the event (16), the following constraint 
is satisfied for all I £ [0,m]: 

(17) 

^max {B{x,ri)) > ri/Tn = Vi. 

We again drop the dependence on n in the notation of u;, and vi is indepen¬ 
dent of n. 

These constraints are simultaneously satisfied if and only if there is a 
strictly increasing sequence of (1 + A:) integers 0 = IqJi, ... ,lk between 0 
and m such that 

Vli+i > 'ymax(B(x,nJ) > if f < A, 

'^max {B{x,ri^)) > Vm- 

We define the following events for all i £ [0,m], with the conventions 
Ik+i = m + l and = 0: 

A-i = {Umax(B(3:,r/J) > 

Bi = > Umax(B(x,D,)) > 

A = {Umax([B(x,D,)] \ [B{x,ri._^)]) > 
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Notice that the events do not depend only on i, but on the whole sequence 
of li- Notation is easier this way. In particular, the former paragraph reads 

m 

(18) '^[L^y> rm\Txy<Tn]< E E 

k=0 

Let us consider the filtration generated by the lines intersecting B(x,r) 
for increasing r, that is = (T(nj 5 (- 3 ,,.)]). Then Ai, Bi and Di are all 

measurable, and Di is independent of J-n. 

Moreover, the difference between Ai and Di is only on the event that 
there is a line faster than hitting B{x,ri^_^). Since k+i — 1 > /i, this 

never happens under So that Ai n = DiCi From this, we 

deduce 


F[Bi\Bo,.. .,Bi_i] < P[^j|.Bo,. ..,Bi_i] since Bi C Ai 

= F[Di\Bo,.. .,Bi_i] since Ai n Bi_i = A n A-i 

(19) 

= P[A] since A indep. of Bj for j <i 
< P[j4j] since Did Ai. 

Recall that the number of lines faster than v hitting a ball of radius r is 
a Poisson variable with parameter with c = /irf([i?(0,1)]). We 

may then compute 

P[^*] = 1 - exp{-crf-\~l]Z^) 

H H+1 ^ 

( 20 ) = 

< since xq < ri. 

_ 2-{7-i)(^i+i-h)p(^p) with 

p(ro) = 2^"^cro"'^A'“\ 

Recalling the convention Ik+i = m + 1, we obtain 

F[Bo n A n • • • n A-i n Ak] 

< P[Ro]IP[A|A] • • ■^[Ak\Bo,Bi,... , A-i] 

k 

<[]P[A] by (19) 

i=0 

k 

< by (20) 
i=0 
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^2-(7-i)(™+i)p(^Q)fc+i. 
Substituting into bound (18), we obtain 

m 

¥[L:,y > rm\T.y <Tn]<J2 E 



k=0 0=lo<-"<l}^<m 

(21) 

= E (T) 

k=0 ^ ' 


= 2-(™+i)(7-i)p(ro)(l+p(ro))”^ 

(22) 

<(2-(7-i)(i+p(^o)))”^+'. 


Let At < 7 — 1. We now choose our free parameters: 

“ (^27-1-k - 1 ) ’ 

so that 

l+p(ro) = 2 ^-l-^ 

m = [(n + 1 )/Atj. 

Substituting into bound (22), we get 

nL^y>r^\T^y<Tr,]<2-^"^+^'>^ 

Since f’[Txy > Tn] < 2 “^”'+^^ we have thus proved that with probability 
at least 1 — 2“"’, the Euclidean length is bounded by 

Lxy ^T’m 

where C is a positive constant depending on T and k, but not on n. 

Thus, Lxy has a (5-moment for all 5 < k. Indeed, 

OO 

HLiy] < Y, 2-”(2(”+i)/^(n + l)V(7-rf)(7)5 

n=0 

OO 

< ^2W'‘-i)”0(n'^/(7-«*)) 
n=0 

< OO if (5 < At. 

Since At is only constrained by At < 7 — 1, we have a 5-moment for all 5 < 7 — 1. 
Since 7 > (i > 2 , this completes the proof. □ 
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Remark 5.1. The main part of the proof is really just saying that the 
Euclidean diameter of a Il-ball has finite mean. Together with Theorem 3.1, 
this implies that the random metric space generated by the R-length is 
almost surely homeomorphic to 

Notice that for k close to 7 — 1, the main term in the sum (21) is the 
one for k = 0. In other words, it is easier to have one single extremely fast 
line close to x than to have many successively faster lines, if we want to be 
abnormally fast at each distance of x. 

A first way to improve on the moments of L^y starts with noticing that if 
we use a line to go far away very fast, we need to use another line to come 
back, since a geodesic never crosses itself. 

More precisely, the author conjectures that the structure yielding long 
geodesics with highest probability is the following: a line with speed at least 
2rjTxy hits R(x,e), another hits B{y,e) and they hit a common e-ball at 
distance r. Since the two first events have probability in and the 

other in we would conclude that K[Lx'y^^~^] = oo, but < oo 

for all 5 < 27 -|- d — 3. 

6. Finite intensity of long-distance network. We now turn to property 4 
of a SIRSN, that is their key property. 

Intuitively, this means that the SIRSN contains “highways.” If we look at 
all the geodesics simultaneously, truncating each geodesic by deleting balls 
around its endpoints, their total length in each compact set is finite: the 
geodesics largely re-use the same segments in each region. Contrast with the 
Euclidean case where the whole space is used. 

In our context, we have to prove that the intensity p(I) of the following 
long-distance network B is finite: let {H„,n € N*} be a collection of Poisson 
processes on with intensity n times Lebesgue, all independent from R, 
and coupled so that C Write H = IJneN* Then 

T'= IJ {gxy\{B{x,l)U B{y,l))). 
x,yeE. 

Notice that almost surely Qxy is unique for all x and y, since the dense 
point set H is countable. 

By translational invariance, it is enough to prove that the intersection 
of J- with a given ball has finite mean Hausdorff measure of dimension 1. 
Indeed, if ^ = mi{F n B{x,r)), then E[^] = —p(I). Notice that scale 

equivariance yields similar results if the long-distance network F is defined 
by removing balls of any fixed radius instead of radius I. 

The main argument relies on the pigeon-hole principle: a geodesic getting 
close to a prescribed point must use one of a very few fast lines close to 
that point, and must use them again to draw away. And by uniqueness of 
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geodesics, two geodesics with two common points must agree between those 
points. 

Theorem 6.1. Let 7 > d > 2. With the above notation, let i = mi{J- n 
B{t), i)) he the length of the long-distance network in B{0, |). For e < Smax; 
with probability at least 1 — e, this length is less than C(ln(C'i/e))^, for 
constants Smax; C and Ci depending only on 7 and d. Consequently, there 
is a finite moment of exponential form: for any 6 < \fC, 

(23) E[exp((5\/Z)] < 00 . 

In particular, i has finite mean. 

Proof. Since H is countable, by Theorem 4.8, almost surely all geodesics 
between its points are unique. In the proof, we only consider such geodesics 
and sub-geodesics, so that we assume uniqueness. 

Consider the balls B{0,^) and i?(0, |). We call their set difference i?(0, |) \ 
B{0, |) the border. 

Now if a point x G i?(0, |), then i?(0, |) C B{x, 1). Hence, geodesics with 
an endpoint in H(0, |) make no contribution to i. We have 

l<mi( IJ (g,„ynB(o,^ 

Hence, geodesics gxy making contributions to i are structured in the fol¬ 
lowing way, illustrated in Figure 4. 

• They hit H(0, |) for the hrst time at point s on the corresponding sphere. 



Fig. 4. The geodesic from x to y hits B{0, |) for the first time at s and the last time at 
z. It hits B{0, |) for the first time at u and the last time at v. 
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• Then they hit B{0, for the hrst time at point u on the corresponding 
sphere. 

• Then they hit for the last time at point v on the corresponding 

sphere. 

• Then they hit for the last time at point z on the corresponding 

sphere. 


In particular, the contribution to i is included in the sub-geodesic guv 
Now, uniformly on x and y, the time Tgz between s and 2 : is bounded 
by the time diameter of B{0, |). So that, by Theorem 3.1, with probability 
1 — e/ 2 , we have 


(24) 


Tsz < Te 


= T, 


2/3,1 


In ■ 


e • Sr, 


i/(7-i) 


We call fast lines the lines faster than = 1/{6T^). Conversely, lines 
slower than are slow lines. We write = •Sve ni?( 0 , |) for the intersection 
of all those fast lines with the ball 5(0, |). 

The number of fast lines hitting i?(0, |) is a Poisson variable with param¬ 
eter Ag = — = C' 2 ln(C'i/e) for explicit constants Ci and C 2 , with 

2 'Ug 

Cl = 2lerasx ^ 2. We recall that the moment generating function of such a 
Poisson variable X is E[e*^] = exp(Ae(e* — 1)), and use Chernoff bound to 
get 

P[A: > CaAe] < 

= exp(A£(C '3 - 1 - Cs InCs)) with e* = Cs 

\ C2(C73-l-C3lnC3) 




with Cs chosen big enough to have C' 2 (C '3 — 1 — Cs InCs) < — 1 . 

Hence, with probability at least 1 — |, there are at most C 3 A£ fast lines 
hitting i?(0, |). With probability at least 1 — e both this event and Tgz < 7/ 
are true. We assume both from now on. 

Since Tgz < 7/, the intersection of gsz with slow lines has length at most 
VelTe = \. 

So that, since s and u (resp., v and z) are at least | apart, the geodesic 
must have fast segments for length at least ^ between s and u (resp., v and 
z), that is, 


mi{Vengsu) > 
mi{Veriguz) > I- 
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Now 


so that 


< #{fast lines} • Diam 


(•(•■I 


< 4 ^ 
“ 3 



mi 0 mi(Ve (8) Vs) 
mi (g) mi {qsu n Vs (8) Qyz n V^) 


<^5 



2 


Hence, by the pigeon-hole property as illustrated in Figure 5, we may find 
a maximal family {g^} of at most (C' 5 (ln^)^ -|- 1 ) geodesics such that: 

• 5 ® is a geodesic between x® and y® in H \ - 6 ( 0 , |), passing through s®, u®, 
u® and z® dehned as for s,u,v and z above. 

• Any geodesic gxy contributing to i crosses one of the 5 ® when hrst and 
last crossing the border, that is, there are points t and w such that t G 
9s^u^ n gsu uud w G g^i^^ hi gvz- 



Fig. 5. Close parallel curves in the figure agree. They are separated to make the figure 
more readable. The three thin black lines are the fast lines. The light dashed lines and 
black dotted curves are a family of geodesics {g'}. The solid curve connecting x and y is 
the geodesic g^y. The geodesic gxy has a common point t with the black dotted geodesic 
on a fast (black) line when first crossing the border, and another w when last crossing the 
border. Hence, they agree between u and v. Any other geodesic contributing to I would 
meet one of the geodesics in the family {g'} in the same way. 







32 


J. KAHN 


By uniqueness of geodesics, 5 * and gxy coincide on gtw In particular, they 
coincide on guv Hence, the intersection TCi B{0, is included in the finite 
number of geodesics . We may then conclude by separating contributions 
from fast and slow lines: 

i <'^mi{guiui) 
i 

< mi ( Ve) + ^ mi [guiui n (H \ Vg)) 

<c(in!h)'. 

Theorem 6.2. The network M made of the time geodesics is a SIRSN. 

Proof. Property 1 of a SIRSN is a consequence of almost sure unique¬ 
ness of geodesics between two points, that is Theorem 4 in dimension at 
least 3, and Theorem 4.4. in Kendall’s (2015) article in dimension 2. 

Property 2 of a SIRSN is because the underlying Poisson line process is 
invariant by translation and rotation. As for change of scale, the underlying 
Poisson line process is invariant by a transformation where scale is multiplied 
by a and speed by Hence, all paths have their time length 

multiplied by the same so that the geodesics are the same and 

M is invariant. 

Property 3 of a SIRSN is Theorem 5.1. 

Property 4 of a SIRSN is Theorem 6.1. □ 

7. Conclusion. We have established that the improper Poisson line pro¬ 
cess with adequate speed limits yield a SIRSN. 

Along the way, a few questions have been raised. Is there an easier, more 
natural way to prove uniqueness of geodesics? What are the tightest mo¬ 
ments of the Euclidean length of a geodesic? When can we generalize this 
construction using geodesics from a random geodesic metric space? 

On a more general note, we may wonder which property of our network 
M translate to general SIRSNs, or to SIRSNs made of geodesics of a metric 
space. For example, it should be easy to show that not all SIRSNs have the 
equivalent of Lemma 4.6 or Corollary 4.9: it is certainly not true for the 
binary hierarchy model by Aldous. 

We might also raise a few typical questions in stochastic geometry. Is there 
only one geodesic connecting a prescribed point x to infinity, like in dimen¬ 
sion two? In many models, infinite geodesics have an asymptotic direction. 
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For SIRSNs, this property looks unlikely, and characterizing the random 
walk of the angle as a function of the Il-distance to x looks worthwhile. 
What is the law of a typical cell in the tessellation generated by the network 
connecting the points of Hi, an intensity 1 Poisson point process? 

Finally, another somewhat tangential direction of research would be to 
study more closely the properties of the random metric space. For exam¬ 
ple, being a SIRSN entails coalescence of geodesics, a very hyperbolic-like 
property. We may also draw comparisons with a well-known random metric 
space such as the Brownian map [see, e.g., Le Gall (2014)]. 

The Brownian map is a random metric space homeomorphic to the sphere 
S^. It has Hausdorff dimension 4. All its geodesics minus their endpoints is of 
Hausdorff dimension 1. The cut-locus of its geodesics starting from a given 
point has Hausdorff dimension 2, and the topology of an open continuous 
tree. 

On the other hand, we have shown that our metric space is homeomorphic 
to It should be easy, by scaling arguments, to show that its Hausdorff 
dimension is (dy — d)/{'y — d), which is bigger than d. In particular, with 
d = 2 and 7 = 3, we have the same dimensions as the Brownian map. If any 
geodesic can be appropriately approximated by geodesics between points of 
Poisson point processes, it should also be easy to show that the geodesics 
minus their endpoints is of dimension 1. However, the cut-locus might have 
a very different behaviour. 

Acknowledgements. I would like to thank Wilfrid Kendall for introduc¬ 
ing us to the problem during a talk, and helpful discussions since then. 
Remarks from a referee have greatly improved the presentation of the pa¬ 
per. 
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